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Abstract — Phonon dispersions in (100) silicon nanowires 
(SiNW) are modeled using a Modified Valence Force Field 
(MVFF) method based on atomistic force constants. The 
model replicates the bulk Si phonon dispersion very well. 
In SiNWs, apart from four acoustic like branches, a lot of 
flat branches appear indicating strong phonon confinement in 
these nanowires and strongly affecting their lattice properties. 
The sound velocity (Vsnd) and the lattice thermal conductance 
(ki) decrease as the wire cross-section size is reduced whereas 
the specific heat (Cv) increases due to increased phonon 
confinement and surface-to-volume ratio (SVR). 

I. Introduction 

Silicon nanowires (SiNWs) are playing a vital role in 
areas ranging from CMOS |Q to thermo-electricity 1 2 1 . The 
finite extent and increased surface-to-volume ratio (SVR) in 
these nanowires result in very different phonon dispersions 
compared to the bulk materials This work investi- 

gates the effect of geometrical confinement on the phonon 
dispersion, the sound velocity (Vsnd), the specific heat (Cy), 
and the lattice thermal conductance (k,i) in (100) SiNWs 
using a Modified Valence Force Field (MVFF) model based 
on atomistic force constants. 

Previous theoretical works have reported the calcula- 
tion of phonon dispersions in SiNWs using a continuum 
elastic model and Boltzmann transport equation ||4j, atom- 
istic first principle methods like DFPT (Density Functional 
Perturbation Theory) |5J-^8J and atomistic frozen phonon 
approaches like Keating- VFF (KVFF) 0, 1^21 • Thermal 
conductivity in SiNWs has been studied previously using 
the KVFF model |11|, |12|. The continuum approaches fail 
to capture the atomistic effects like anisotropy in Cy for 
nanowires, etc., whereas the first principle models cannot 
be extended to very large structures due to their heavy 
computational requirement. The VFF model overcomes 
these shortcoming by capturing the atomistic effects and 
solving phonon dispersions and transport in realistic struc- 
tures 1^, |[TT|. However, a simple KVFF model does not 
reproduce the bulk dispersions very well and hence cannot 
be used confidently for nanostructures |13|. The MVFF 
model overcomes the verifiable shortcomings in bulk. We 
use the MVFF model for phonon calculations in SiNWs 
and believe them to be more robust and predictive than the 
KVFF model |13|. 

This paper is organized in the following sections. Section 
|Il| briefly presents the phonon model and its application to 
SiNWs. Also the computation of different wire lattice prop- 
erties is outlined there. The results on phonon dispersion and 
other lattice properties are presented in Sec. [Ill] Conclusions 
are summarized in Sec. ITvl 




Fig. 1. The short range interactions used for Phonon dispersion calcu- 
lations. The interactions are (a) bond-stretching(Q;), (b) bond-bending(^) 
(c) bond stretch coupling(7) (d) bond bending- stretching coupling (6) (e) 
coplanar bond bending coupling(A). 

II. Theory and Approach 

A. Modified VFF method 

The Valence Force Field (VFF) model is a force constant 
based atomic potential (U) calculation method. The MVFF 
model we present here is based on the combination of two 
different extended VFF models, (a) VFF model from Sui 
et. al. 1 14] which is suited for non-polar materials like Si, 
Ge and (b) VFF model from Zunger et. al. p5| which is 
suited for polar materials like GaP, GaAs, etc. These meth- 
ods have successfully modelled the phonon dispersion in 
various semiconductor materials |[T4|-|[T6|. The force con- 
stants represent the various kinds of interaction between the 
atoms. For phonon modeling in diamond (and zinc-blende) 
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Fig. 2. Projected unitcell structure of a (100) oriented Silicon nanowire. 
The white atoms (for visual guidance) show the freely vibrating surface 
atoms. Grey colored atoms are inside the wire 
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Fig. 3. Comparison of the MVFF method (lines) and the experimental 
data (dots) ^17J for bulk Silicon phonon dispersion at 80K. 



lattices we consider the following interactions (Fig. [T]), (a) 
bond-stretching (a), (b) bond-bending (c) bond stretch 
coupling (7) (d) bond bending- stretching coupling (S) and 
(e) co-planar bond bending coupling (A). The first two terms 
(a and b in Fig. [T]) are from the original Keating model 
|[16 |, which are not sufficient to reproduce the bulk phonon 
dispersion accurately in zinc-blende semiconductors fTS] , 
|Pl-| . Hence, higher order interactions (term c, d and e 
in Fig. [T]) are included in the original KVFF model. The 
force constants in the MVFF model depend on the bond 
length and angle variations from the ideal value |14|. Due 
to this reason MVFF is also known as the ' quasi- anharmonic 
model'. 

The motion of atoms in the semiconductor structure 
is captured by a dynamical matrix (DM). The dynamical 
matrix is assembled using the second derivative (Hessian) 
of the crystal potential energy (U). The bulk DM has 
periodic boundary conditions along all the directions (Born- 
von Karman condition, Fig. [T]) due to the assumed infinite 
material extent along all the directions. For free-standing 
SiNW, periodic boundary conditions are applied only along 
the length of the wire (X-axis) assuming an infinitely long 
wire, whereas the surface atoms (Y-Z axis) are free to 
vibrate as illustrated in Fig. [2]). The complete calculation 
details are provided in Ref. 1 13 |. 

B. Lattice property calculations 

A wealth of information can be extracted from the phonon 
spectrum of solids. The description of the few important 
ones are provided below. 

Sound Velocity: One important parameter is the group 
velocity (Vgrp) of the acoustic branches of the phonon 
dispersion. Near the Brillouin Zone (BZ) center Vgrp charac- 
terizes the velocity of sound (Vsnd) in the solid. Depending 
on the acoustic phonon branch used for the calculation 
of Vgrp, the sound velocity can be either (a) longitudinal 
(Ysnd^i) or (b) transverse (Ysnd.t)- Thus, Vsnd is given by. 
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where, A and q are the phonon polarization (for longitudinal 
or transverse direction) and wave vector respectively. 

Lattice thermal properties: In a semiconductor where its 
two ends are maintained at a small temperature difference 
(AT), its specific heat (Cy) and thermal conductance (ki) 



can be evaluated from the phonon dispersion using the Lan- 
dauer's formula |[T8| . The temperature dependent specific 
heat is given by |T9|, 
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where ks, K T, and n are Boltzmann's constant, reduced 
Planck's constant, mean temperature, and number of sub- 
bands, respectively. The temperature dependent lattice ther- 
mal conductance (hii) for a ID conductor is given by |j6J, 
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where, M(cj) is the number of modes at a given frequency 
uj. In the next section we present the results on phonon 
dispersion in SiNWs and the lattice thermal properties. 

III. Results and Discussion 

A. Experimental benchmark 

The MVFF model has been bench marked by calculating 
the phonon dispersions in different bulk zinc-blende semi- 
conductors. One of the representative calculation is shown 
in Fig. [3] for bulk Si. The model accurately reproduces the 
important features of the experimental phonon dispersion of 
Si (at 80K) 1 17] in the entire BZ. The force constant values 
for Si are listed in Table |l| The bond stretch coupling pa- 
rameter (7) is zero for Si. This term contributes significantly 
for splitting the longitudinal optical (LO) and longitudinal 
acoustic (LA) phonon branches near the BZ edge only in 
III-V materials like GaAs, GaP |[T5|. 

B. Nanowire phonon dispersions 

The phonon dispersion is calculated in free-standing 
(100) SiNW where the surface atoms are allowed to vibrate 
freely (Fig. |2]). The dispersion in SiNWs with cross-section 
size (W=H) (a) 2nm, (b) 4nm, and (c) 6nm are compared 
in Fig. |4] The number of sub-bands increases as the wire 
cross-section size increases due to an increased number of 
atoms (Fig.|4]). This results in a higher degree of vibrational 
freedom. A comparison of the number of sub-bands within 
10 meV of energy for the different SiNWs considered here 
shows that the 2nm x 2nm structure has ^ 40 sub-bands, 
the 4nm x 4nm has ^ 90 sub-bands and the 6nm x 6nm 
has nearly 110 sub-bands (Fig.|4]). The increased separation 
of phonon sub-bands in smaller wires is a consequence of 
the geometrical confinement. The higher sub-bands have 
mixed acoustic and optical branch like properties. Many of 
the sub-bands in this regime have group velocities close to 
zero which again indicates a strong phonon confinement in 
thin nanowires. In SiNWs these flat phonon branches appear 



TABLE I 

Force constants used for Silicon in MVFF model 
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Fig. 4. Effect of cross-section size on the phonon dispersion in (100) 
SiNW. The wire dimensions are (W = H) (a) 2nm (b) 4nm, and (c) 6nm. 
As the cross-section size increases the number of phonon modes within a 
given energy range increase as shown here. The geometrical quantization 
increases the separation of bands in smaller wires. 

due to the zone-folding of acoustic branches from the bulk 
phonon-dispersion (Fig. [3]). 

Eigen modes of oscillation: The first four vibrational 
Eigen modes are shown in Fig. |5] (arrows show the displace- 
ment direction). The first two modes (Fig. [5] a, b) are the 
flexural modes responsible for NW bending perpendicular 
to the NW axis. The Eigen frequencies for the lowest two 
phonon bands near the zone center are cx q'^. The next 
mode (Fig. |5] c) is the longitudinal mode. The mode no. 
4 is torsional mode which represents the extra degree of 
rotational movement in free standing wires (Fig.|5]d). Mode 
no. 3 and 4 behave like the bulk acoustic branches (oc q near 
the zone center) [IQJ . Mode 1 and 2 appear due to the free- 
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Fig. 6. Sound Velocity (Ysnd) in (100> SiNW with free boun dary . 
As a reference the bulk Vsnd is shown along the (100) direction 1^0]. 
Increasing the wire diameter increases the Vsnd since the vibrational 
degree of freedom increases. 

Standing nature of the nanowires which are not available in 
bulk. Hence, new modes as well as new phonon dispersion 
branches appear in SiNW which are very different from their 
bulk counterpart. This also has a strong impact on the lattice 
properties of the wires which is discussed in the following 
part. 

C. Lattice properties in SiNW 

Sound velocity variation (Vsnd)-' The sound velocity 
(Vsnd) is another way to compare the vibrational modes in 
SiNWs. Vsnd is obtained using the sub-bands 3 and 4 (trans- 
verse and longitudinal modes) of the phonon dispersion. 
Vsnd shows a reduction with decreasing wire cross-section 
size (acoustic mode softening) fTO| (Fig. |6]). This happens 
since the acoustic branches get flatter in small wires due to 
phonon confinement. This strongly affects the thermal prop- 
erties of the SiNWs since the acoustic branches are mainly 
responsible for heat conduction. A smaller sound velocity 
also results in smaller thermal conductance (discussed later). 

SiNW specific heat (Cy) variation: The specific heat 
of silicon nanowires increases with their decreasing cross- 
section size (Fig. [7]). This can be attributed to two reasons, 
(i) phonon confinement due to small cross-section and (ii) 
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Fig. 5. Eigen modes of the atomic displacement shown in a SiNW with 
W, H = 1.63nm. The atomic displacement are shown at ~ for (a) 
mode = 1, (b) mode = 2. These modes are cdXXtd flexural modes and cause 
bending of the NW. (c) mode = 3, called the longitudinal mode causes 
atomic motion along the wire axis, (d) mode = 4, called the torsional 
mode, observed only in wires, represents the extra rotational degree of 
freedom. 



O 



720 



700 



680 



660 



640 




Temp = 300K 



^'bulk = ^^S"* @ ^^^^ 



3 4 5 

W = H [nm] 



Fig. 7. Variation of the specific heat {Cy) in Silicon nanowires with 
width. As a reference the specific heat for bulk Si is shown by a dashed 
line. SiNWs have larger Cy than bulk Si |7 |. The specific heat increases 
as the cross-section size of the wire decreases. 
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Fig. 8. Variation in the ballistic lattice thermal conductance (ki) (a) 
with temperature for different cross-section size SiNW and (b) at room 
temperature (300K) for different NW width size. 

an increased surface-to- volume ratio (SVR) in smaller wires. 
Geometrical confinement separates the phonon bands in en- 
ergy with decreasing wire cross-section size (Fig. |4]) which 
makes only the few lower energy bands active at a given 
temperature (see Eq[2]). This explains the increasing energy 
needed to raise the temperature of the smaller wires. Also 
the increasing SVR results in a higher partial phonon DOS 
associated with the wire surface, which further enhances the 
specific heat with decreasing wire cross-section |7|. 

Lattice thermal conductance (k,i) variation: The 
modified phonon dispersion in SiNWs significantly affects 
the ballistic hzi |11|. The ballistic hzi increases with 
temperature as shown in Fig. [8^. The spread of the 
Bose-Einstein distribution (Fbd = [exp(ficj/A;^T) — 
in energy increases with increasing temperature. This 
results in a larger number of modes contributing to the 
heat transfer ||6| (see Eq. ([3])). This explains the increased 
ballistic hzi at elevated temperatures. Larger wires show 
more hii since, (i) they have more phonon sub-bands 
resulting in more number of modes (see Eq. ([3])) and (ii) 
they have a higher acoustic velocity which is responsible 
for larger heat conduction (Fig. |6]). The ballistic hzi of 2nm 
X 2nm Si nanowire reduces by a factor of 8 as compared 
to the 6nm x 6nm SiNW (Fig. 

An important point to note here is that k,i is expected to 
decrease further in smaller wires due to phonon scattering 
by other phonons, interfaces and boundaries |11| which is 
neglected in the present study. The main message here is 
that even ballistic phonons show a reduction in thermal 
conductance with narrow SiNWs, which comes from (i) a 
modification of the phonon dispersion (ii) phonon confine- 
ment effects in SiNWs. 

IV. Conclusions 

A modified atomistic force constant based model has been 
developed for the calculation of phonon dispersion in zinc- 
blende lattices. This model shows a very good agreement 
with experimental bulk phonon data. An important impact of 
device miniaturization is the presence of surfaces (and hence 
higher surface-to-volume ratio) which affects the phonon 
dispersion in SiNW. The clear demarcation of phonon 
bands as acoustic and optical branches become vague in 
SiNW. This manifests itself in the form of modified lattice 
properties in SiNWs compared to bulk Si. We find that 
decreasing wire cross-section reduces the sound velocity and 



the lattice thermal conductance whereas the specific heat 
increases. The geometry dependent calculation of phonon 
dispersion in small nanowires is therefore important to 
properly understand their thermal properties. 
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